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Flavour-singlet pseudoscalar and scalar mesons 

UKQCD Collaboration, C. Michael, M. S. Foster* and C. McNeile 

Theoretical Physics Division, Dept. of Mathematical Sciences, University of Liverpool, Liverpool L69 
3BX, UK 

We measure correlations appropriate for favour-singlet mesons using dynamical quark configurations from 
UKQCD. Improved methods of evaluating the disconnected quark diagrams are presented. The scalar and pseu- 
doscalar meson channels are explored, with special reference to the mixing of scalar mesons with scalar glueballs. 



1. INTRODUCTION 

To study flavour singlet mesons, we need to 
consider quark loops which are disconnected (of- 
ten called hairpins), namely evaluate TrFM"^ 
where the sum (trace) is over all space (for zero 
momentum) at a given time value and all colours 
and spins. Here M is the lattice fermion matrix, 
r is a combination of the appropriate 7-matrix 
and a product of gauge links if a non-local oper- 
ator is used for the meson. 

Using a random volume source ^ (where 
(C*?) = 1 for the same colour, spin and space- 
time component and zero otherwise) then solving 
Mcj) = ^, one can evaluate unbiased estimates of 
the propagator M~y from {i*cj)x) where the av- 
erage is over the Ns samples of the stochastic 
source. The drawback of this approach is that 
the variance on these estimates can be very large, 
so that typically hundreds of samples are needed. 
Here we present a method which succeeds in re- 
ducing this variance substantially at rather small 
computational expense. 

The variance reduction is based on expressing 
the fcrmion matrix M as 

M^C + D = C{1 + C-^D) = (1 + DC-^)C{1) 

where C is easy to invert, for example the SW- 
clover term which is local in space. Then we have 
the exact identity 

Af^i = C-^ - C-^DC'-^ + ... + {-C~^D)"'C-^ 
+ {-C-^D)"^M~\-DC-^)"' (2) 

*present address; MVC, Manchester Computing, Manch- 
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with ni+n2 = n = m + l. 

Using the stochastic estimate for on the 

rhs will reduce the variance of the estimate given 
by the Ihs since the terms not involving can 
be evaluated either exactly (for example terms 
with odd powers of D vanish in the evaluation 
of a local trace) or as a subsidiary stochastic cal- 
culation with more samples since no inversion is 
required. 

A special case of this (ni = 712 = 2) with Wil- 
son fermions (for which C = 1 and the terms with 
up to 3 powers of D vanish for TrM"'^) was used 
by the bermion group previously. 

Using larger values of ni and 712 implies that 
the estimate of AI~^ is very non-local. To eval- 
uate correlators between traces at ti and ^2, one 
must require that the samples of stochastic vol- 
ume source used in the two cases are different so 
that there is no bias. We use Ns = 24 stochas- 
tic samples and this condition is readily imple- 
mented. This number of samples was chosen to 
make the stochastic sampling error smaller (actu- 
ally about 50% in the cases to be discussed next) 
than the intrinsic variance from one gauge config- 
uration to another. The computational effort is 
equivalent to that in obtaining two conventional 
propagators (from two sources of all colour-spins) 

We now apply this technique to 12'^ 24 lattices 
with Nf = 2 at P = 5.2 with Csw = 1-76 from 
UKQCD H]. We use valence quark masses equal 
to the sea-quark mass determined by the k- values 
as given in Table 1. In order to improve the statis- 
tics we measure the disconnected diagrams on 
configurations separated by less trajectories than 
for the connected correlators giving 253, 169 con- 
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figurations at k = 0.1395, 0.1398, respectively. 

2. PSEUDOSCALAR MESONS 

Since we are exploring Nf ~ 2 we describe the 
flavour non-singlet (7 = 1) pseudoscalar meson 
as TT and the flavour singlet (/ = 0) as i]. It is 
possible to apply our results to the physical case 
of 3 light quark flavours and the rj, rj' splitting as 
has been done from the quenched case |^ , but we 
do not discuss this here. 

We use local and fuzzed operators matching 
the calculation of the connected meson propaga- 
tors ||]. Then including our results for the hair- 
pin diagrams with the appropriate sign, we can fit 
the flavour singlet correlator to the rj mass and ex- 
cited states. Even though, as discussed above, we 
have measured the disconnected quark loop cor- 
relator from every time slice to every other with 
a stochastic error which is effectively negligible, 
the signal is very noisy. This arises since the dis- 
connected part of the r] correlation has an error 
which does not decrease with increasing time sep- 
aration T, unlike the connected correlator where 
the error is roughly constant as a percentage of 
the signal with increasing t. To illustrate this, we 
show the ratio of disconnected D to connected 
component C of the 77 correlator in flg. 1. 

The line shows the expectation coming from 
one state (tt) contributing to C and one state 
(77) contributing to C — D with a mass difference 
in lattice units of 0.16. There is some sign in 
our data of the impact of unitarity for matched 
valence and sea quarks: namely that D/C ap- 
proaches 1 from below at large t. Making a fit 
with this constraint to D/C, it is possible to esti- 
mate the r] — TT mass difference and we see that it 
does increase as the sea quark mass is decreased. 
Assuming that the rj mass is constant as the quark 
mass is decreased, we would get an r] mass in the 
chiral limit of around 800 MeV, of course with an 
uncontrolled systematic error 

Because the flavour singlet correlator becomes 
noisy so rapidly with increasing time separation, 
it is very valuable to use a more extended basis 
for the meson (eg. local and fuzzed operators) to 
help in extracting the ground state signal. The 2 
state fit from t = 3 to 9 is shown in Table 1. 
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Table 1 



The zero-momentum correlator has contribu- 
tions increasing like for spatial volume 
whereas the noise increases like L^. Thus we find 
a better signal to noise ratio for smaller volumes - 
hence we present here only the results for L = 12. 
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Figure 1. Disconnected to connected ratio. 

3. SCALAR MESONS 

For scalar mesons, we have the interesting sit- 
uation that the qq mesons and the glueball can 
mix. Within the quenched approximation, it is 
possible to estimate this mixing Q|. We make 
a preliminary study in the quenched approxima- 
tion at /3 = 5.7 of 100 configurations on 12^^ 
24 lattices with SW-clover valence quarks having 
K = 0.14077 with Csw = 1-57 (here a « 1 GeV 
and the quark mass is close to strange). We use 
4 scalar meson operators (i) closed Wilson loops 
(glueball operators) of two different sizes (Teper- 
smeared) and (ii) qq operators which are local 
and separated by fuzzed links. Using the glueball 
mass of 0.97(4) from ref. and scalar qq mass of 
1.48(15) from fitting the connected correlations, 
we are able to fit the disconnected correlation and 
the hairpin- Wilson loop correlation from i-values 
of 1 to 3 with a mixing given by Ea = 0.4 as- 
suming Nf = 2. This fit (see fig. 2) assumes only 
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Table 2 

ground state contributions so that the systematic 
error on E from this assumption is hard to esti- 
mate. The mixing estimated by ref Q is similar 
in magnitude at these lattice parameters {E « 0.3 
GeV) but they claim that on extrapolation to the 
continuum limit a much smaller value is obtained. 
Here we are using clover improvement so order a 
effects are suppressed. If our quenched mixing 
strength were be to applied to the scalar mass 
matrix, it results in a downward shift for N 
of the lattice glueball mass by 20%. 



1.5 



0^^ /S = 5.7 12-'24 Np=0 
7t/p = 0.65 Ero=1.2 

G: Wilson loop 
D: Quark loop 

LrGr 

HH: Connected quark 




0.0 



Figure 2. Quenched scalar correlations. 

We can now measure directly the scalar spec- 
trum for Nf = 2 to explore this. We again use 
both glueball and qq operators and fit with 2 
states the 4x4 matrix of correlations for t from 
2 to 10: the results is shown as mpso- in Table 
2. We find that the mass obtained from fitting 
only the glueball correlations (moB) is consistent 
with the full fit, as it should be. Moreover, we see 
a surprisingly low scalar mass - as emphasised in 
fig. 3 which compares with quenched results and 
the SESAM Nf^2 values ||. 

We do expect a relatively light flavour-singlet 
scalar mass because of mixing effects as described 
above which would reduce the mass by 20%. This 
could explain in part our low scalar mass but 



other explanations are also worth exploring. For 
example the order corrections might be anoma- 
lously large for our lattice implementation (e.g. 
twice as large as in the quenched Wilson case). 

Another possible explanation of the light 
flavour-singlet scalar mass we find (comparable 
to our pion mass) would be a partial restoration 
of chiral symmetry in a finite volume. Our spatial 
size is 1.7 fm and no evidence of finite size effects 
was seen in a study of flavour non-singlet corre- 
lators We have made a preliminary study of 
flavour singlet correlators on 16'^ spatial lattices 
to check for flnite size effects and, although the 
signal from the larger spatial volume is relatively 
noisier, we do see some evidence (at the 2a level) 
of a higher scalar mass on the larger spatial lat- 
tice. 
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Figure 3. The scalar mass versus a^. 
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